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Abstract 

The pp — > ppir differential cross section has been measured with the ANKE spec- 
trometer at COSY-Jiilich for pion cms angles between 0° and 15.4° at a proton 
beam energy of 0.8 GeV. The selection of diproton pairs with an excitation energy 
E pp < 3 MeV ensures that the final pp system is dominantly in the spin-singlet 
x 5o state. The kinematics are therefore very similar to those of pp — > di: + but with 
different spin and isospin transitions. The results will thus provide a crucial extra 
test of pion production models in nucleon-nucleon collisions. 

The cross sections, which are over two orders of magnitude smaller than those 
of pp — > dir + , show a forward dip, even stronger than that seen at lower energies. 
This behaviour is well reproduced in a theoretical model that includes P-wave AN 
states. 
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Single pion production in nucleon-nucleon collisions, NN — > NNir, is the first 
inelastic process that can be used to test our understanding of the underlying 
meson-baryon dynamics of the NN interaction [1,2,3]. By far the cleanest re- 
action to study is pp — > dn + , where the differential cross section and multitude 
of spin observables that have been measured over the years [4] can confront 
the different theoretical models. 

In contrast, very little was known about the pp — > ppir reaction, though the 
unexpectedly large it production rate observed near threshold [5] led to a 
flurry of theoretical activity [6]. Now in cases where the excitation energy E pp 
of the final protons is very small, due to the Pauli principle, this reaction will 
excite only the J p = + (^o) diproton state. Despite having kinematics very 
similar to those of pp — > dii + , this reaction involves different transitions in the 
NN system and, in particular, the role of the A isobar is expected to be much 
suppressed because the S-wave AiV intermediate state is forbidden. 

Some information on the transitions involved has been extracted from quasi- 
free pion absorption on 3 He [7]. However, previous measurements of the dif- 
ferential cross sections for pp — > {pp} s n° have only been carried out up to 
a beam energy of 425 MeV, with a cut imposed on the excitation energy of 
Epp < 3 MeV [8,9]. This value satisfies the requirement that the spin-singlet 
S'-wave (^o) final state, here denoted by {pp} s , should dominate while pro- 
viding reasonable statistics. One important feature of the experimental data 
is that with the E pp selection the cross sections show a forward dip whereas, if 
no cut is applied on the excitation energy, then for the higher beam energies 
there is a forward maximum [9,10]. 

Now 425 MeV is below the threshold of NA production; the data show no sign 
of being influenced in any clear way by the A and we need to go to higher en- 
ergy to investigate the effects of P-wave AiV systems. As part of a programme 
to investigate the small E pp region in intermediate energy nuclear reactions, in 
particular in large momentum transfer deuteron breakup reactions [11,12], we 
have carried out a high statistics measurement of the pp — > {pp} s n° reaction 
at T p = 800 MeV for pion cm angles below 15.4°. 

The experiment was performed at the magnetic spectrometer ANKE [13], 
placed at an internal target position of the COSY COoler SYnchrotron [14]. 
Fast charged particles, resulting from the interaction of the proton beam with 
the hydrogen cluster-jet target [15], were registered in the Forward Detector 
(FD) system [16]. Its hodoscope provided a trigger signal and an energy-loss 
measurement. It also allowed a determination of the differences in arrival times 
for particle pairs hitting different counter elements. The tracking system gave 
a momentum resolution a p /p ~ 0.8 - 1.2% for protons in the range (0.5 - 
1.2) GeV/c. 
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Fig. 1. Scatter plot of the magnitudes of the momenta of two charged particles 
detected in the FD. The selection procedure introduces a slight bias as to which 
particle is called "1", but this does not affect the subsequent analysis. 

The trigger used required the crossing of the two planes of the scintillation 
hodoscopes by at least one charged particle but, in the subsequent off-line 
analysis of the pp — > {pp} s vr° reaction, only events with two tracks in the FD 
were retained. In Fig. 1 is shown a two-dimensional scatter plot of the mag- 
nitudes of their two momenta corresponding to about half of our statistics. 
Due to the limited angular acceptance of ANKE, there are kinematic correla- 
tions for reactions with two and three particle in the final state. One therefore 
sees in the figure islands corresponding to pp — > dn + and bands resulting from 
pp — > pnn + and pp — > ppn°. Candidates for the latter reaction are well sepa- 
rated from the other processes. Furthermore, in approximately 75% of cases 
the particles hit different counters in the hodoscope and the difference in their 
arrival time could also be used in the selection. The dir + pairs coming from the 
pp — > dn + reaction, which could potentially provide the most serious physical 
background, are separated from the pp pairs from pp — > ppn in time difference 
by more than at 8 ns, whereas the actual resolution is better than 0.5 ns 

The distributions of missing mass squared, M|, are shown separately in Fig. 2 
for single-counter and double-counter candidates with low excitation energy in 
the pp system, E pp < 3MeV. In both cases one sees a very clean tt° peak centred 
at 0.021 (GeV/c 2 ) 2 , which agrees with m 2 to well within our experimental 
precision. The widths of the Gaussian fits are compatible with those obtained 
from Monte Carlo simulations; the marginally narrower peak in the single- 
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Fig. 2. Distributions in the square of the missing mass for candidates for the 
pp — > ppX reaction with excitation energy E pp < 3 MeV and 9^ < 15.4° when 
the protons (a) hit different counters, and (b) the same counter. From the indicated 
positions of the 7r° peak and the 2ir° threshold it is seen that single and double pion 
production can be clearly separated. Gaussian fits to the 7r° peak plus a constant 
background yield a total number of 7r° events in (a) and (b) of respectively 4425 
and 1008. 

counter data is due to these events generally having a smaller opening angle 
resulting in the kinematics being slightly better denned. The backgrounds are 
small and slowly varying and two-pion production can be clearly excluded in 
either case. There is a small excess of events observed on the left side of the n° 
peak. These may correspond to single photon production through pp — > {pp} s 7 
and so the regions indicated by dashed lines in Fig. 2 have not been included 
in the Gaussian fits. Since this interpretation is not unambiguous and these 
events might still correspond to good 7r° events, we have added an extra 2% to 
the systematic error. Given that the two data sets are completely compatible, 
they have been grouped together in the subsequent analysis. The resolution 
in excitation energy for the combined pp — > ppn° events was cr(E pp ) 0.2- 
0.3 MeV for E pp in the range 0-3 MeV. 

The value of the luminosity needed to determine the cross section was found 
by comparing the yield of pp elastic scattering, measured simultaneously with 
the other reactions, with that deduced from the SAID data base [17]. The 
integrated luminosity obtained from this is L- mt = (6.72 ± 0.26) x 10 34 cm -2 , 
where the error comes mainly from averaging over the angular bins. 

In order to determine the triply differential cross section d 3 a / (dili™ dE pp ), 
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Fig. 3. (a) Number of pp — > pp7r° events in the interval 8.9° < < 13° as a function 
of E pp ; (b) The same data corrected for acceptance and detection efficiency and 
presented as differential cross sections. Only statistical errors are shown. The curve 
results from passing the Migdal-Watson function |T(.E PP )| 2 of Eq. (1), multiplied 
by phase space, through a Monte Carlo simulation of the ANKE apparatus and 
normalising the predictions to the summed experimental histogram. Similar results 
are found for the other angular intervals. 

events selected in the range < E pp < 3 MeV were divided into groups of 
equal intervals in cos 9^. The energy spectrum of counts in the angular in- 
terval 8.9° < 9p™ < 13° is shown in Fig. 3a. The spectra in the other intervals 
demonstrate a similar behaviour. Values of the corresponding cross sections 
were obtained from such a distribution by taking into account geometrical 
acceptance, efficiency in the track recognition algorithm for two particles, in- 
teractions in the constituent materials, efficiency and resolution of the detec- 
tors and other known effects on the basis of a Monte Carlo simulation of the 
ANKE setup. This leads to the histogram with the statistical errors presented 
in Fig. 3b. 

The rapid rise of the spectrum with E pp from threshold illustrated in Fig. 3b 
is typical of all intermediate energy reactions where one produces proton pairs 
and is induced by the pp final state interaction. We have indeed observed ex- 
actly the same phenomenon in the pd — * (pp)n reaction with the same appara- 
tus at ANKE [11,12]. This effect is often parameterised, in the Migdal-Watson 
approximation [18], by the square of the low energy pp elastic scattering am- 
plitude for which 
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where |C(7y)| is the Coulomb penetration factor evaluated at rj = am p /2k 
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Fig. 4. Distribution of acceptance-corrected pp — > 7r° events with < 3MeV 
over cos#*, where 9* is the angle between the directions of the proton and diproton 
momenta in the centre of mass of the diproton. Note that the vertical scale does 
not start from zero. 



a^m p /E P p/2, and 5 is the combined Coulomb-nuclear phase shift. This has 
been evaluated numerically for the Reid soft core potential, for which the 
scattering length a pp = —7.8 fm. 

The Migdal-Watson factor of Eq. (1) was used as an event generator together 
with phase space to provide candidates which were then traced through the 
experimental setup, taking into account all its known features. The resulting 
smoothed curve, shown in Fig. 3a, provides a semi-quantitative description 
of the data which is quite sufficient for our purpose, where we quote cross 
sections summed over energy. The data are a little above the curves at the 
higher E pp and we cannot exclude some small P-wave contribution though 
globally the angular distribution of the pp system in its rest frame shown in 
Fig. 4 is consistent with isotropy. It should be noted that the 1 P final state 
would also produce a flat distribution. 

Due to the identity of the initial protons, the differential cross section is an 
even function of cos#£ m and in Fig. 5 it is plotted versus cos 2 9^ m . The results 
show a monotonic decrease towards the forward direction and, as seen from the 
figure, they can be well parameterised by the linear function a(l + b sin 2 #£ m ), 
where a = (704 ± 22 stat ± 32 syst ) nb/sr and b = 5.6 ± 1.2. With the same E pp 
cut as used here, a similar forward dip was observed in this reaction at lower 
energies, T p < 425 MeV [8,9], though for these energies it was found that b 
was much smaller, being always less than 1.4. Since, for such small values of 
E pp , the final diproton must be dominantly in an S-wave, constraints from 
spin-parity and Fermi statistics then require the pion to be in an even partial 
wave. As a consequence, the forward dip was attributed to an interference 
between the pion s and d waves [9]. Given that the influence of rf-waves might 
be expected to increase with energy, it is perhaps not surprising that we find 
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Fig. 5. The measured pp — > {pp} s ir° differential cross section for E pp < 3 MeV as a 
function of cos 2 (#£ m ). The curve is a straight-line fit to the data. 

a larger slope parameter at 800 MeV. 

Preliminary theoretical predictions have been made for the pp — > {pp} s 7r° dif- 
ferential cross section at 800 MeV in a model that includes contributions from 
P-wave AN intermediate states [19,20]. The overall magnitude is similar to 
that which we have observed and, in particular, the forward slope, driven by 
the pion rf-wave, is well reproduced. It is expected that our data, combined 
with those at lower energies, will allow such models to be refined. 

Though we have argued that the kinematics of pp — > {pp} s n° and pp — > dn + 
are quite similar, the underlying dynamics must be very different. This is 
illustrated in Table 1, where we show the values of the two differential cross 
sections and their ratio R(tt° / 7r + ) in the forward direction obtained at different 
energies. The results seem to indicate that there might be a broad minimum 
in R in the A region of the pp — > d-K + reaction. 

The GEM collaboration has recently published high precision data on the 
ratio of the forward production of pions in the pp — ^ ir + d and pp — > 7i + pn 
reactions at 981 MeV [21]. Since the spin-singlet pn final state interaction 
has a much sharper energy dependence than that of the triplet, from the 
shape of the pion momentum spectrum they could put an upper limit on the 
amount of pn singlet produced. Integrating over excitation energies E pn < 

3 MeV, it is seen that the ratio of -r^;(pp — > {pn} s ir + ) to 7777 (j>£> — > dir + ) could 

be at most about 5%. If Coulomb effects are ignored, the cross sections for 
spin-singlet production through pp — > {pn} s n + and pp — > {pp} s n° should be 
identical. Though Coulomb suppression in the {pp} s final state will be large, it 
looks very doubtful whether the study of the spectrum alone will be sufficient 
to isolate the pp — > {pn} s n + cross section in view of the values presented in 
Table 1. Measuring the proton and pion in coincidence, as has been done for 
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Table 1 

Zero degree differential cross sections for pp — > {pp} s ir° with < 3 MeV from 
the present experiment at 800 MeV with the lower energy data being taken from 
Ref. [9]. The values of the pp — > dir + cross sections are obtained from the SAID 
SP96 solution with the range of the other solutions being taken as a rough estimate 
of the error bars [17]. The ratio R(tt°/tt + ) of the two pion-production cross sections 
is also presented. 

example in Refs. [22,23] and analysed in a model-independent way [24], still 
only provides upper bounds. The study of 7r° production therefore seems to 
be the most realistic way of investigating the final state here. 

Data on quasi-free pion production in the pd — > {pp} s X reaction were ob- 
tained as a by-product of our deuteron break-up measurements [11] and 
these results will be interpreted in terms of the sum of the cross sections 
for pp — > {pp} s 7T° and pn — > {pp} s 'K~ ■ At 800 MeV it will then be possible 
to subtract the 7r° contribution reported here in order to obtain data on 7r~ 
production. 

It is intriguing to note that a very similar ratio to that of Table 1 has been ob- 
served for backward dinucleon production in the pd — > {pp} s n and pd — > dp re- 
actions at intermediate energies [11]. Now such a connection would be natural 
within a one-pion-exchange mechanism, where the large momentum transfer 
pd — > dp reaction is driven by a pp — > dn + sub-process [25,26]. More quanti- 
tative estimates of the pd — > {pp} s n cross section, where the pp — > {pp} s 7r° 
sub-process is used rather than the pp — > dn + , are currently under way [27]. 

It is seen from Table 1 that there is a real lack of data on the pp — > {pp} s vr° re- 
action in the A region and this could be usefully filled by further experiments 
at ANKE. It should also be noted that, unlike the complicated spin struc- 
ture connected with the pp — > dn + reaction, only two spin amplitudes which 
are functions of cos 2 9 n are required to describe the pp — > {pp} s n° reaction. 
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These can be isolated, up to an unmeasurable overall phase, by determining 
the proton analysing power and the initial pp spin correlation C xx . Both of 
these experiments can be carried out at small angles using ANKE [28] and 
the resulting amplitude analysis will tie down even further ttNN dynamics at 
intermediate energies. 

This work was supported in part by the BMBF grants ANKE COSY-JINR, 
Kaz-02/001 and Heisenberg-Landau programme. We are grateful to many 
other members of the ANKE collaboration who provided strong support for the 
measurement. Important discussions with J. A. Niskanen and correspondence 
with I. Strakovsky are also gratefully acknowledged. 
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